Is D+j(2632) the first radial excitation of D*(2112)? 
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Abstract 

We present a quantitative analysis of the D^(2632) observed by SELEX mainly focusing on the 
assumption that D,tj(2632) is the first radial excitation of the 1~ ground state Z)*(2112). By solving 
the instantaneous Bethe-Salpeter equation, we obtain the mass 2658 ± 15 MeV for the first excited 
state, which is about 26 MeV heavier than the experimental value 2632 ± 1.7 MeV. By means of 
PC AC and low-energy theorem we calculate the transition matrix elements and obtain the decay 
widths: T(D+j -» D+n) = 4.07 ± 0.34 MeV, T(D+j -» D°K + ) ~ T(r(L>+ J -> D+K°) = 8.9 ± 1.2 
MeV, and the ratio T(Djj -> D°K + )/T(D+j -> D+rf) = 2.2 ± 0.2 as well. This ratio is quite 
different from the SELEX data 0.14 ± 0.06. The summed decay width of those three channels 
is approximately 21.7 MeV, already larger than the observed bound for the full width (< 17 
MeV). Furthermore, assuming 1)^(2632) is 1" state, we also explore the possibility of S — D wave 
mixing to explain the SELEX observation. Based on our analysis, we suspect that it is too early to 
conclude that Z)jj(2632) is the first radial excitation of the 1~ ground state D*(2112). More precise 
measurements of the relative ratios and the total decay width are urgently required especially for 
S — D wave mixing. 
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I. INTRODUCTION 



Last year the SELEX Collaboration reported the first observation of a narrow charmed 
meson D+^2632) with the mass of 2632.5 ± 1.7 MeV and decay width T tot < 17 MeV at 90% 
confidence level This state has been seen in two decay modes, D^rj and D°K + , with 
a relative branching ratio D°K + )/T(D+j D+rj) = 0.14 ± 0.06. This relative 

branching ratio is rather unusual because the phase space of channel D°K + is 1.5 times 
larger than that of the channel D^rj and according to the SU(3) flavor symmetry the ratio 
T(D+j -»• D°K + )/T(D+j -> D+rj) > 2.3 2] if the quark content of ££,(2632) is cs, and 
the reference |3J gave an estimate of T(D^j — > D°K + )/T(D^j — > D+rj) = 7.0. Since the 
/},t 7 (2632) is above the threshold for D°K + and Dfr], and its decays to these final states are 
Okubo-Zweig-Iizuka (OZI) allowed processes, it is also strange that it has a narrow width 
(< 17 MeV). 

Why this state is so narrow and dominated by the Dtp decay mode has inspired a lot 
of theoretical interest 

BBflflBBBBQQQQflQ 

, with many possible 

explanations, e.g. the first radial excitation of -D*(2112), a tetraquark, a hybrid, a two- 
meson molecule, a diquark-antidiquark bound state, etc. Among them the most plausible 
candidate is the first radial excitation of .D*(2112), the 2 3 5*i state, since it is the likely 
conventional meson with component of cs which has not been found. And even though the 
mass of 2632 MeV is about 100 MeV below the traditional potential model prediction (for 
example, see Q]), the similar situation has been found in the new narrow charm-strange + 
meson D s j(2317), favored as the conventional cs state. Furthermore, Simonov and Tjon Q] 
showed a possible mechanism to shift down the meson mass by the coupled-channel analysis 
method. Chao P gave a reasonable argument that the node structure in the 2 3 5*i state may 
explain the small decay width and unusual decay modes. 

In this letter, following the references , we assume _Dj}(2632) as the first radial 
excitation of vector meson D*(2112). Then we solve the instantaneous Bethe-Salpeter equa- 
tion |l7j], and give our prediction of mass for 2 3 S' 1 state. At the same time, we obtain 
the wave functions of the relevant mesons. By using the reduction formula, PCAC rela- 
tion and low energy theorem, we write the transition S-matrix as a formula involving the 
light meson decay constant and the corresponding axial current transition matrix element 
between two heavy mesons, which in turn can be written as an overlapping integral of the 
relevant wave functions. We give a detailed consideration on the node structure in 2 3 S' 1 wave 
function. We further try to make a rough estimate on the full decay width and the ratio 
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T(Djj -> D°K + )/T(Djj -> D+77) by assuming D+j(2632) as the D-wave 1" vector state as 
well as considering the possible S — D wave mixing. 

II. PREDICTION OF MASS 

In our previous papers we solved exactly the instantaneous Bethe-Salpeter equa- 

tions (or the Salpeter equations 20j) to describe the behavior of CT and 1~ states. In Refs. 
1^ | . we have given the mass spectra and wave functions of pseudoscalar D s , D° and D + , 
and the mass spectra of all these heavy 0~ states, which are consistent with experimental 
data quite well. For the instantaneous Bethe-Salpeter equation of 1~ states [19|, we have 
used the same parameters as used in [18:], only changing the parameter Vq to fit data be- 
cause Vq is the only parameter which can be different for 0~ and 1~ states when their quark 
contents are the same. In this letter we first solve the full Salpeter 1~ equation through 
fitting the ground state mass at 2112 MeV (£)*(2112)) to fix the parameter Vq. After fixing 
the parameter Vq, we then give the predictions of the first radial excitation state. Since we 
are only interested in the leading order estimate, which lets us see the problem clearly, we 
assume in our calculation that the first radial excitation state Dfj is a pure S wave state, 
excluding the possible mixing between S wave and D wave states. Later in this letter we 
will also consider the D wave contribution. Our prediction for the mass of the first radial 
excitation state is 2658 ± 15 MeV, which is about 26 MeV heavier than the experimental 
data 2632 ± 1.7 MeV. 



III. TRANSITION MATRIX ELEMENT 



In this section we give a detailed description for the calculation of the transition matrix 
element. By using the reduction formula, the transition S'-matrix of decay Dfj — ► D°K + 
(similar to the channel Dfj — > D + K°) can be written as: 

(D(P n )K(P f2 )\Djj(P)} = J d*xe iP f- x (Ml - Pj 2 )(D(P n )\$ K (x)\Dtj(P)), (1) 

where P, Pfi and Pf 2 are the total momenta of initial state 2 3 S\, final state heavy meson 
and final state light meson, respectively (see Figure 1); M# is the mass of the final light 
meson. The PCAC relates the meson field <&k( x ) with a current <&k{x) = M i f <9 M (g7^7ss), 
where f K is the decay constant of light meson; q is u for K + , d for K°, respectively. Then 
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FIG. 1: Diagram for the transition of -Dj}(2632) decay, 
the S'-matrix becomes 



{D(P fl )K(P f2 )\Dfj(P)) 



(Ml - P%) 



M 2 K f R 



J d 4 xe iP f^(D(P fl )\d^q W )\Dfj(P)) 



M 2 K f K 



J d 4 xe iP f-*(D(P fl )\q W \Dfj(P)), 



and with the low energy theorem this equation can be written approximately as: 



-^E f d'xe^%D{P }1 )\q W s\Dtj{P t )). 

J K J 



(2) 



(3) 



And, finally, we have 



(DiP^KiP^Dfjim = (2nmP t - P fl - P f2 



— ?'P M 
—Jl, 

f K 



D(P fl )\q lfll5 s\Dtj(m. (4) 



For the decay Dfj — > D+r), the equation is similar, but there is rj — rj' mixing. Including 
the rj — rj' mixing, the PCAC relation reads as 

® v (x) = cos9$ m (x) + sm9® vo (x) 
- C ° s6> gM( ^75W + ^75^-2s7 M 75S ^ | sing ^ ^7^75^ + ^75^ + ^755 ^ 



= ^(57^758) 



rjo-IVO 



-2 cos e 



+ 



sing 



(5) 



where f m , f Vo are the decay constants of octet r/s and singlet r/ , respectively. Then the 
transition matrix element becomes 

(DsiPnHPf^Dijim = (2 7 r)V(P i - P fl - P f2 )(-iP? 2 ) 

4 



X 



-2M„ 2 cos# M 2 sin# 

1 _|_ 1 



(6) 



Note that the decays D+j(2632) -> D+ are OZI-rule-allowed, so we may ignore anomaly of 
the PCAC safely in Eqs. (Effijl . 

According to the Mandelstam formalis m |2ll| . at the leading order, the matrix element 
(D(Pfi)\q'y fl 'y 5 s\D^j(Pi)) can be written as [22]: 

dq 



(D(P fl )\q W \Dtj(P)) 



(27T) 



-Tr 



— — tO-tt^C^TmTs 

?™c + m a M * 



(7) 



where m c and m q is the mass of c and g (q = u, d, s) quark; r is the recoil three dimensional 
momentum of the final state D meson. p c , p s are the momenta of c and s quark in the 
initial meson, respectively. And the relative momentum q is defined as q = p c — m m ^ m P 

1 //, . in the center of mass system of the initial meson, and it becomes q = p c = —p s ; 



m c +m a 

y?+ + , and tp 



n 



are called the positive energy wave function of the initial and final heavy 



mesons, and ip~ 
in paper jl8| |: 



n 



— 7o(v5^ + ) + 7o- The positive energy wave function of D° can be found 



M 



n 



■/I 



uj c + uj q q*(m c - m q 

; ^To- 
me + m a 



+ 



m c + m q \ 



m q u! c + m c u) q (m q uj c + m c uj q ) 



(8) 



where oj c = \jvn? c + q 2, and u q = yTn q + g 2 ; <fi(q), f2(o) are the radial part wave functions, 
and their numerical values can be obtained by solving the full Salpter equation of 0~ state. 

Since the wave function <f ++ {q) directly relates to the first radially excited state, there 
is only one node in it. We note that the node structure can possibly play an important 
role to explain both the narrow decay width and the dominance of Dfr] decay mode, to 
a certain degree. The value of the radial wave function of 2 3 Si state becomes negative as 
a function of relative momentum q when it crosses the node (denoted as go), that is when 
\q\ > go- Therefore, this negative part of the wave function tends to cancel the contribution 
to the decay width from the positive part of the same wave function, which can be the 
reason to have a small decay width. In the overlapping integral of the wave functions Eq. 



-r, and this will even strengthen the negative 



(171) there is the shift of momentum, q 

contribution because the shift will move the peak of the radial wave function of the final 
D meson close to the node. We note that the value of the overlapping integral will be 
significantly affected by the momentum shift as the larger momentum shift will cause the 
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larger negative contribution, and in turn the smaller overlapping integral (smaller decay 
width). Since the recoil momentum \f\ of D meson is larger than that of D s , about 1.43 
times, and the momentum shift in case of DK is 1.56 times larger than that of the D s r\ 
channel, larger negative contribution results in DK channel than in the D s 7] channel, which 
can be the very reason for the small ratio of the branching ratios DK/D s rj. 

As already pointed out in Ref. j3j, there may be the S — D mixing in 1~ state , so 

we will try to estimate the possible mixing contributions here too. We first assume that this 
1~ state is a pure S wave state 2 3 Si, and later we will consider the case of D wave state. 
The wave function of 1~ can be written as jlo| : 

ip Pi (q) = M { + V>27o) 7o [709 ■ *l>i + M~ ?• e)lfe] ^ + f " ^ 



(m c + m q )q 



(u c + uj q )q 2 J 



and the positive energy wave function is defined as: 



^(O) = A tlo<P Pi (<?)7oA+ = T^i^clo + m c + #)7o¥^ (q)lo^(^ g lo -m q -fl, (10) 

where e is the polarization vector of a 3 Si state; A+, A+ are the energy projection operators 
for quark and antiquark; the numerical value of the radial part wave function ipi(q), ip2{q) 
will be obtained by solving the exact Salpter equation of 1~ state. 

IV. NUMERICAL RESULTS AND DISCUSSIONS 



For our numerical calculations we use the following 



parameters: m c = 1755.3 MeV, 
is! U = 130.7 MeV f m = 1.26/ w 
24( | . The mass of octet rjs can be 



m s = 487 MeV, m d = 311 MeV, m u = 305 MeV 
MeV, f Vo = 1.07 U MeV |23j and f K = 159.8 MeV 
estimated by the Gell-Mann-Okubo relation 3m 2 g = 4m 2 K — m 2 , the mass of singlet rjo can 
be obtained by m 2 o = m 2 +m 2 , — m 2 g , and the mixing angle can be estimated by the relation 
m 2 o = sin 2 6*m 2 + cos 2 6*m 2 ,. [So our input values are m m = 564.3 MeV, m Vo = 948.1 MeV, 
e = -9.95°.] 

We estimate the mass of the first radial excitation of D*(2112) by solving the full Salpeter 
equation of l - state, and our prediction is 2658 ± 15 MeV, which is 26 MeV heavier than 
the central value 2632 ±1.7 MeV measured by SELEX. The theoretical uncertainties are 
estimated by varying all the input parameters simultaneously within ±5%, and we choose 
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the largest possible error. We obtained the wave functions of the corresponding states, and 
we show radial part wave functions <fi, y?2 of D s in Figure 2; <pi, ip 2 of D° in Figure 3; ipi, 
ip2 of 2 3 Si in Figure 4. We show the overlapping part of 2 3 S'i wave function i^i(\q\) and D s 
wave function <pi(\q\ — \r n l) 5 when it has the biggest momentum shift — ^ — |r n I i n 

r J- vi i| m c +m s I D a \n oo m c +m s I D s \ 

Figure 5. We also show i^i{\q\) and D° wave function (fii(\q\ — m m + m V D S) i n Figure 6. One 
can see that in 2?S\ wave function the node localizes at go — 0.54 MeV, and the biggest 
momentum shifts are mc \r" n \ = 0.447 MeV for D°, and mc |r" I = 0.283 MeV for Dt. 

rric+mu 1 d° 1 ' rn c +'m s 1 D s 1 s 

The decay width for the two body final state can be written as: 

87rM 2 ' 1 1 ; 

Here the matrix element T for D+j -> D°K + and D+j -> D + K° is 

T = ^-(DCP/OI^Tfeap+^PO), (12) 

and for — > it is 



T = P% 



-2M%cos9 M*sm9 



(DiPrifavWlD+W). (13) 



We obtain the following decay widths: 

T(D+j -> D+77) = 4.07 ± 0.34 MeV, (14) 

T(Djj -> 73°^+) = 8.87 ± 1.23 MeV, (15) 
T(D+ J -> JJ + ^°) = 8.74 ± 1.22 MeV. (16) 
And the corresponding ratios are 

r(D+j -> J D°^ + )/r(/j+ - jj+7?) ~ rp+ - D + K°yr(Dtj - ~ 2.2 ± 0.2. (17) 

As can be seen, we cannot find the dominance of D s r), and the ratio is much different from the 
experimental ratio 0.14 ± 0.06. The sum of those three decay widths is 21.7 MeV, already 
larger than the observed bound for full width (< 17 MeV) given by SELEX experiment. 
Therefore, our assumption of the pure S wave vector state cannot fit the experimental data. 

Before considering the S — D mixing in 1~ vector state, let us first assume _Djj(2632) to 
be a pure D wave state. In this case, we do not solve the full Salpeter equation, but instead, 
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FIG. 2: Radial wave function y?i(|*1) and v?2(|*1) of Df. 
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FIG. 3: Radial wave function </?i(|*1) and </?2(|*1) of D°. 
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FIG. 4: Radial wave function V'id'Zl) an d ^2 (1*1) of 2 3 5i. 
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FIG. 5: The overlapping part of radial wave function and fi(\q\ — m ™ c m I ? Ds I ) when the 

momentum of tpi has the biggest shift — ^ — \r n I. 
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FIG. 6: The overlapping part of radial wave function V'iG'Tl) an d Vid^l 
momentum of ip\ has the biggest shift 



m c +m u 1 D° 1 



when the 



m c +m u I d° ' ' 



we give the leading order estimate for the pure D wave explanation. The wave function of 
a pure D wave state (S — L coupling) can be simplified as: 



^(9) = J2 M ( x + io)A s z)i>nu z {.q) < ^s z -,ii z 1 lj z > , 



(18) 



where 1 = 2. With this wave function as input and by using the same parameters as for the 
S wave state, we can solve the approximate Salpeter equation, i.e., only the positive energy 
part solution is considered. The mass of this pure D wave state is 2672 MeV, and the radial 
wave function of this pure D wave state is shown in Figure 7. 

To estimate the matrix element in this case, we calculate the overlapping integral of initial 
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FIG. 7: Radial wave function ip(\q\) of £>+ J (2632) if it is a I? wave state. 



and final wave functions. There is only one term, which can give a non-zero contribution: 



J [ZTrr m r , + m n 



(19) 



m c + m q 

where i — 1, 2; c is a constant; /(g 2 ) is a shorthand notation for the terms which are functions 
of (f; e al3 (l z ) is the polarization tensor, which couples with the polarization vector e 7 (s 2 ) to 
make the total polarization vector e a (j z ): 



£ e a (s z )e^(l z ) < ls z ;ll z | lj z >= 



o p/S p-y 



pa 07 pa p/3 



■ (« = 2) 



(20) 



M 2 7 VJ " /V * M 2 

Then, we obtain the numerical values for the decay widths when Dfj(2632) is a pure D wave 
state (J p = I - ): 

r(D+j -> D+77) = 3.88 ± 0.32 MeV, (21) 
r(D+ J -> D°fs: + ) = 23.3 ± 3.2 MeV, (22) 
T(D+j -> D + K°) = 21.5 ± 3.1 MeV. (23) 
And the corresponding ratios are 

rp+j - L> x + )/rp+ - l>+77) ~ r(L>+ - L> + x°)/r(L>+ - ~ e.o ± 0.5 , (24) 



which show that the pure D wave assumption does not fit the experimental data obviously. 
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FIG. 8: The decay width with different mixing angle 9, when initial state mass is 2632 MeV and 
it is an S — D mixing state. 



We now consider the possibility of the S — D mixing with a mixing angle 9: 

\D* aJ ) = cos 9\2S) + sin9\lD} , 
\D* s *j) = -sm9\2S) + cos0|lD) . 



(25) 



D°K + and D+j -> D + K° 



We show the decay widths for the channels D^j — > Djrj, D^j 
as a function of the mixing angle 9 in Figure |H1 Numerically it is easy to find that when the 
mixing angle is taken to be 9 = 0.80 77T, the relative branching ratio 



T(Dfj -> D°K + )/r(Dfj -> Df V ) = 0.14 



can fit the SELEX data. We also predict other relative branching ratios 



(26) 



and 



T(D+ - D+K°)/T(D+ - Din) ^ 0.05, 
F(D+ - D+K°)/r(D+ -> D°K+) ~ 0.37, 



(27) 



with a rather small total decay width 



T = r(£>+ -> + r(L>+ -> D ^) + r(L>+ -> D+K ) + • • • ~ 0.35 MeV, (28) 

which is quite smaller than the experimental bound for the full width (< 17 MeV). It is inter- 
esting to note that in fact there is another solution for the mixing angle, 9 = 0.8377T, which 
can also fit the data of the relative branching ratio. For this solution, accordingly, the relative 
branching ratios F(D+j -> D + K°)/V(D+j -> D+rj) ~ 0.20 and F(D+j -> D + K°)/T(D+j -> 
D°K + ) ~ 1.4 are predicted with a slightly larger total decay width T ~ 0.82 MeV. The two 
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solutions present a different relative aspect on the two decay channels Dfj — > D°K + and 

As a final note, we only estimated the OZI allowed strong decay widths D^j — > D s r] 
and D^j — > DK assuming that Dfj is an 1~ vector state. However, there are some other 
mechanisms through which the 1~ vector state can decay to the same final states, which 
we do not consider here. For example, there is the OZI forbidden mechanism for Dfj 
decaying to D s r) through the anomaly term in PCAC, i.e. it can decay through two gluons, 
Dfj — > D s + gg, gg — > rj, which is a decay similar to ^' — > J/^rj (with decay width about 
8.9 keV). As is well known, the OZI forbidden processes are small and cannot improve 
sufficiently the relative ratio to agree with the SELEX data. Therefore, we do not include 
the processes in our analysis. 

As summary, according to the data from SELEX and our estimates made here, it is not 
likely that the _Djj(2632) is the pure S"-wave 1~ vector state. Only the S — D wave mixing as- 
sumption may explain the relative branching ratio T(D^j — > D°K + )/T(D^j — > Dfrj) = 0.14, 
whereas the assumption seems to give a rather small total decay width (to be compared 
with the experimental bound T < 17 MeV). Generally speaking, our results (based on 
PCAC and low energy theorem) are qualitatively similar to those of Ref. |3( (based on 3 P 
model). We find that the node structure for the first radially excited state of the 1~ ground 
state D*(2112) can have a favorable effect on the relative ratio of the decay widths but it 
can not be so big as indicated by the data. The S — D wave mixing may help to explain 
the data, however, more precise measurements on the relative ratio T ^'*i* D „+ . - (and/or 

r(DJj^Dlri) y i 

+ ^ + ) as well as the total decay width are crucial. Therefore, the explanation of 

v sj > 

.Djj(2632) as the first radially excited state of the 1~ ground state D*(2112) is too early to 
conclude. 
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